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I. INTRODUCTION 

During the past decades the Casimir effect [1] found 
a large number of applications in different branches 
of physics (see monograph Q and references therein). 
Among them the applications should be mentioned in 
statistical physics, in elementary particle physics (e.g. in 
the bag model in QCD or in Kaluza-Klein theories) and 
in the cosmology of the early Universe. From the point 
of view of quantum field theory the Casimir effect is a 
specific type of vacuum polarization which appears by 
quantizing the theory in restricted volumes or in spaces 
with non-trivial topology. This polarization results from 
a change of the spectrum of zero-point oscillations in the 
presence of nontrivial boundary conditions. For the case 
of the electromagnetic vacuum between two uncharged 
metallic boundaries, separated by a small gap a, the 
Casimir effect leads to the appearance of an attractive 
force acting on them depending on a and on the funda- 
mental constants h and c only. Such attractive force may 
be alternatively explained as a retarded van der Waals 
force between the two bodies whose conducting surfaces 
are responsible for the Casimir force. 

The Casimir force was firstly measured by Sparnaay 
H for metallic surfaces. For dielectric bodies the corre- 
sponding forces had been measured more frequently, see 
[4-7]. The relative error in the force measurements was 
about 100% in § and in the range of (10-20)% in [4-7]. 

As it was shown in Refs. [8,9] the Casimir force be- 
tween macroscopic bodies is very sensitive to the pres- 
ence of additional hypothetical interactions predicted by 
unified gauge theories, supersymmetry and supergravity. 
According to these theories interactions of two atoms 
arise due to the exchange of light or massless elemen- 
tary particles between them (for example, scalar axions, 
graviphotons, dilatons, arions and others). Their effec- 
tive interatomic potentials may be described by Yukawa- 
and power-laws. After the integration over the volume 
of two macro-bodies one obtains more complicated laws 
for their interaction potentials. It was rather unexpected 
that quite strong constraints for the characteristic con- 
stants of these laws may be found from the experiments 
on Casimir force measurements. 

The constraints under consideration may be found also 
from other precision experiments, e.g., from Eotvos-, 
Galileo- and Cavendish- type experiments, from the mea- 
surements of the van der Waals forces, transition proba- 
bilities in exotic atoms etc (for a collection of references 
on long-range hypothetical forces see JlO|] ) . 

According to the results of [8,11,12] the Casimir force 
measurements of [4-7] lead to the strongest constraints 
on the constants of Yukawa-type hypothetical interac- 
tions with a range of action of 10 -8 m < A < 10~ 4 m. For 
A < 10~ 8 m the best constraints follow from the measure- 
ments of van der Waals forces in atomic force microscopy 
and of transition probabilities in exotic atoms |13| . For 
A > 10 _4 m, as it was shown in [11,12], they follow from 



Cavendish- and Eotvos-type experiments [14-17]. 

In (9) the constraints were obtained from the Casimir 
effect on the constants of power-law potentials decreas- 
ing with distance as r~ n . For n = 2, 3 and 4 they turn 
out to be the best ones up to 1987 (compare (l8)). In 
|]l9f a bit stronger constraints on the power-law interac- 
tion were obtained from the Cavendish-type experiment 
ofRef. @. 

Recently, a new experiment was performed p0| on the 
measurement of the Casimir force between two metal- 
lized surfaces of a flat disc and a spherical lens. The 
absolute error of the force measurements in ||(| was 
AF 10~ n 7V for distances a between the disc and the 
lens in the range 1 /im < a < 6 /im. This was the first 
measurement of the Casimir force between metallic sur- 
faces after || . For the distance a ~ 1 /jm the value of 
the Casimir force in the configuration under considera- 
tion is Fc ~ —3.1 x 10~ 10 N. This means that the relative 
error of the force measurement at 1 /im in [pof is about 
S ~ 3% (note that with increasing a the value of S in- 
creases quickly). In |2^] the active surfaces of the disc and 
the lens were covered by thin layers of copper and gold. 
The use of heavier metals for the test bodies is preferable 
for obtaining stronger constraints on hypothetical inter- 
actions. This follows from the fact that the value of the 
hypothetical forces increases proportionally to the square 
of the density. 

The aim of the present paper is to give an accurate 
calculation of different hypothetical forces which might 
appear in the configuration used in the experiment P0| . 
Also, the corrections to the Casimir force due to distor- 
tions of the surfaces, to edge effects, finite conductivity 
and non-zero temperature will be analysed. On this base 
new constraints on the hypothetical interactions will be 
reliably calculated, which follow from the results of [ p0| , 
as well as their possible improvement. The correspond- 
ing constraints which result for the masses of light ele- 
mentary particles are also discussed. Some preliminary 
results of this kind were obtained in ||l| for Yukawa-type 
interactions and in ]22| for power-law forces. However, in 
[21,22] the corrections to the Casimir force were not dis- 
cussed and different possibilities suggested by experiment 
[20] were not accounted for in full detail. 

The paper is organized as follows. In Sec. II we dis- 
cuss the expression for the Casimir force in a configura- 
tion of a plane disc and a lens and different corrections 
to it are considered taking into account the finiteness of 
the diameter of the disc. Sec. Ill is devoted to the cal- 
culation of Yukawa-type forces in this configuration. In 
Sec. IV analogous results are obtained for the case of 
power-law hypothetical forces. Sec. V contains a care- 
ful derivation of the new constraints on the constants of 
Yukawa- and power-law interactions which follow from 
the experiment performed in p0|| . The possible improve- 
ment of the experiment [^0| is considered in Sec. VI. For 
Yukawa-type interactions it gives the possibility to make 
the constraints about several thousand times stronger. 
This considerable strengthening of constraints may be 
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achieved also for the power-law hypothetical interactions. 
Sec. VII contains the conclusions and some discussions. 
In the Appendix the reader will find a number of mathe- 
matical details concerning the calculation of the Casimir 
and hypothetical forces. 



II. THE CASIMIR FORCE BETWEEN 
A DISC AND A LENS INCLUDING 
CORRECTIONS 

The scheme of the configuration used in the recent 
demonstration of the Casimir force [Q is shown in Fig. 1. 
The Casimir force was measured between the metallized 
surfaces of a flat disc (with radius L = 1.27 cm and thick- 
ness D = 0.5 cm) and a spherical lens (with curvature 
radius R = 11.3 cm and height H = 0.18 cm). The sep- 
aration between them was varied from a = 0.6 fim up to 
6 /im. Both bodies were made out of quartz and covered 
by a continuous layer of copper with A = 0.5 /im thick- 
ness. The surfaces which faces each other were covered 
additionally with a layer of gold of the same thickness. 
Note that the penetration depth of the electromagnetic 
field into gold is approximately <5o ~ 0.08 /im <C A. By 
this reason when calculating the Casimir force one may 
consider the interacting bodies as being made of gold as 
a whole. 

The experimental data obtained in pp|] has been con- 
fronted with the theoretical result 
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This formula is valid for the configuration of a small lens 
situated near the center of a large (strictly speaking, in- 
finite) disc at zero temperature. It was first derived in 
and reobtained by different methods afterwards (see, 
e.g., [2, 23]). In Ref. @, (@) was derived from the well 
known result for two infinite plane parallel plates using 
the Proximity Force Theorem H]. Note, that according 
to our notations the attractive forces are negative and 
repulsive ones are positive. 

Actually, in the experiment j20| the diameter of the 
disc was not much larger, but even smaller than the size 
of the lens. Therefore it is of great interest to calculate 
corrections to Eq. ([!]) due to the finite disc size. For 
this purpose we use the approximation method devel- 
oped earlier for the calculation of the Casimir force and 
which may be applied to the case of two bodies with arbi- 
trarily shaped surfaces [2, 23, 25, 26]. According to this 
method the potential of the Casimir force can be obtained 
by summation of the retarded van der Waals interatomic 
potentials over all pairs of atoms in the bodies under con- 
sideration with a subsequent multiplicative renormaliza- 
tion (the latter takes into account a large amount of the 
non-additivity of these forces). The method was tested 
and successfully applied in the above cited papers. As a 
result the Casimir force may be calculated according to: 
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where r\2 is the distance between the atoms belonging to 
the first and to the second body, i &(ei,£2) is a tabulated 
function depending on the static dielectric permittivities 
of the test bodies (for its explicit form see [2, 26, 27]). 
In the limit of perfectly conducting surfaces (which is 
of interest here) £1,62 ^ 00 and ^(£1,62) — * irhc/24. 
In |^| the relative error of the values given by (g) was 
examined. It was shown to be of the order of 0, 01% 
for configurations which do not much differ from that of 
two plane parallel plates. This is just the case in the 
experiment [^0| , because only the top of the lens and its 
vicinity contribute essentially to the Casimir force. 

The integration in (g) for the configuration of a lens 
and a disc (Fig. 1) may be performed analytically. For 
this reason we introduce a spherical system of coordinates 
with the origin in the curvature center of the lens. The 
angle 8 is counted from the horizontal axis directed out 
of the origin to the left (see Fig. 1). Then Eq. (||) for Uc 
takes the form 
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where the integration limits are defined as follows (prime 
is used for the lens) 
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Using the potential energy in the form of (||) with the 
integration limits (|J) we slightly change the experimental 
configuration converting the disc into a part of a trun- 
cated cone (the corresponding addition to the disc vol- 
ume is AV/V « D/R ps 4%). This increasing takes 
place, however, near the outer boundary of the disc and 
practically does not influence the result for the Casimir 
force. 

The integrals in (3) may be calculated along the lines 
presented in Appendix. Putting k — 3 in (All) we ob- 
tain the result for the power-law interaction with a power 
equal to 2k + 1 = 7: 
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where = min(#i, # 2 ). 

We rewrite (5) in a more convenient form by introduc- 
ing the new variables t = cos'd, x\ — tr\, x 2 — tri and 
using Eq. (A6) for rir 2 : 
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where t = max {R/y/R 2 + L 2 , (R - H)/R). It is seen 
that the value of to depends on the r elative sizes of the 
lens and of the disc. If L < \J2RH one has to w 1 — 
L 2 /(2R 2 ). But if L > V2RH, then t = (R - H)/R and 
the Casimir force does not depend on a further increase 
of L due to the quick decreasing of the retarded van der 
Waals potential with distance. 

Calculating the Casimir force by the first equation of 
(2) the integration with respect to X2 is removed resulting 
in the expression 
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This result coincides with (1) if we consider the limit 
£,£1,2 - ► 00 (infinite disc made of a perfect metal) keep- 
ing the lowest order in the small parameter a/R. Inte- 
grating eq. (7) explicitly we find that the main contri- 
bution to the result, depending on the size of the disc, 
appears at the third order in a/R [28] 
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For the parameters of experiment [20] the inequality 



L < \/2RH holds. In this case it follows from (8) 
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It is easily seen from (|]) that the correction to the 
Casimir force due to the finite size of the disc does not 
exceed its maximal value 6 x 10 - 7 which is achieved 
for a = 6 fim. That is why one actually may neglect 
this correction when confronting the measurements of the 
Casimir force with the theory. 

Let us discuss the corrections to the Casimir force (|Tj) 
which are significant for a small spatial separation of 
the lens and the disc, a w 1 /im. It is reasonable to 
start with the corrections due to the finite conductivity 



of the metal covering the test bodies [2,29-31]. It is well 
known that for a in the micrometer range the penetration 
depth of the electromagnetic field into the metal is fre- 
quency independent and inversely proportional to the ef- 
fective plasma frequency of the electrons: So = c/uj p . For 
two plane parallel metallic plates the corrections for the 
Casimir force due to the finite conductivity were found 
in [29-31] (see also gj). Up to the first two orders in the 
relative penetration depth So /a the result is: 
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Using the Force Proximity Theorem [24] it is not diffi- 
cult to modify (10) for the configuration of a lens and a 
disc: 
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Note that the first order correction in (|ll]) was found 
firstly in |20]| . The plus sign in front of it in [^o| is a 
misprint (this is also clear from general considerations 
according to which Fc(a) is constant in sign for all So 
and tends to zero in the formal limit 5o — ► 00 so that the 
correction should be negative Q). 

For gold, as it was mentioned above, So ~ 0.08 fim and 
for a = 1 /im the correction to the Casimir force F^ due 

to the finite conductivity achieves 23% of F^ \ The be- 
haviour of this correction (in relative units) with increas- 
ing a is shown in Fig. 2 (curve 1). The experimental data 
in p0| do not support the presence of corrections of the 
result ([!]) being so large (let us remind that the relative 
error of the force measurements at 1 /Ltm was about 3% 
and with such an accuracy Eq. (|j) was confirmed). 

According to [M the reason of this contradiction is 
the inapplicability of ([h]), ( |TT| ) for gold (in more detail, 
the approximation for the effective dielectric permittivity 
which was used in [fS0(] to derive the first order correction 



in (10) does not take into account the large imaginary 
part of the refraction index for gold |20|). On the other 
hand the corrections due to the finite conductivity (which 
are in agreement with |30j) were found in |3l]] (see also 
Q) in a more general impedance approach without use 
of the dielectric permittivity. Consequently, Eqs. ( |l0| ) 
and ([ll]) are still valid for gold covered surfaces and the 
contradiction with the experimental data is still present. 
The most reasonable way to resolve this problem is to 
take into account the corrections of the Casimir force due 
to the deviations of the surfaces from the perfect shape 
(in |2(J such corrections were not discussed). 

Let us now assume that the surfaces of the lens and of 
the disc are covered by some distortions with character- 
istic amplitudes A\ and A2, correspondingly. Then the 
general result for the Casimir force up to second order in 
the relative amplitudes of the distortions takes the form 
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where the explicit expressions for the coefficients in terms 
of the functions describing the shape of distortions were 
found in [23]. 

For small stochastic distortions, instead of (12), one 
has [23] 



F c (a) = F^(a)+A d F^(a) 
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where £1,2 are the dispersions of the stochastic perturba- 
tions on the surfaces. The same result is valid for the non- 
stochastic short-scale distortions regardless of the specific 
shape of the functions describing them. Here 6%^ should 
be substituted by A\_ij\pl. Note that the result (13) can 
be obtained using the Proximity Force Theorem ]24[ from 
the corresponding result for two plane parallel plates |32] . 

Speculating that the characteristic values of the disper- 
sion is of order Si « 0.1 pan we get from ( |l3|) a positive 
correction to the Casimir force which is equal to 12% of 
Fq (a) for a = 1 /im. In Fig. 2 (curve 2) the dependence 
of this correction on a is shown in relative units. 

Note that the first order correction resulting from Jl2] ) 
appears only in the case when there are some large-scale 
deviations of boundary surfaces from the perfect shape 
for which the non-perturbed force is calculated. Using a 
realistic estimation of Aj/a w 0.1 for a w 1/im the first 
order correction in (^) may amount as much as 30% 
of F^ [^3|. In the experiment |2(J the radius of the 
lens R = (11.3 ± 0.1) cm was determined with a rather 
large absolute error. Consequently, a large-scale devia- 
tion of the surface of the lens from the perfect shape could 
have taken place. A special inspection of the lens used in 
[po}| . which was not done, is required to determine the ac- 
tual size (and shape) of the large-scale deviations. After 
that it would be possible to calculate the coefficients Cm 
and C01 and the corresponding correction to the Casimir 
force. 

Together with the correction AdF^\a) due to the 
short-scale distortions it may easily compensate the neg- 
ative correction to the Casimir force due to the finite 
conductivity of gold. This is possibly the reason why in 
pof neither corrections to the finite conductivity nor to 
the surface distortions were observed at a » 1 /im. 

One more correction to the Casimir force (Q) is due to 
the non-zero temperature T. It is calculated, e.g., in |3^] 
for two plane parallel metallic plates (see also j| ) . For a 
lens and a disc the corresponding result may be obtained 
by the use of the Proximity Force Theorem and has the 
form 
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Here £ = ksTa/ '(he), kB is Boltzmann's constant and 
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The relative value of the temperature corrections, cal- 
culated with (Jl4|), (15) is shown in Fig. 2 (curve 3) in 
dependence on a. It is seen that for a — 1 /im it is approx- 
imately 2.7% of F^ \ But for a = 6 /im the temperature 
correction is A T i^ 0) = 1.74F£ 0) . 

Fig. 2 will be used in Sees. V, VI which are devoted 
to obtain stronger constraints on the constants of hy- 
pothetical interactions from the Casimir force measure- 
ments (the other forces which may contribute in the ex- 
periment, i.e. electric one should be subtracted to get 
the result for the Casimir force with possible corrections 
discussed above, see Cfl). 



III. THE YUKAWA- TYPE HYPOTHETICAL 
INTERACTION BETWEEN A DISC 
AND A LENS 

The Yukawa potential between two atoms of the inter- 
acting bodies may be represented in the form 



V Y u = -aN 1 N 2 hc — e" 
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-ri 2 /A 
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where a is a dimensionless interaction constant, A = 
h/(mc) is the Compton wavelength of a hypothetical par- 
ticle which is responsible for the rise of new interactions, 
and Ni are the numbers of nucleons in the atomic nuclei; 
they are introduced in (16) to make a independent on 
the sort of atoms. 

The potential energy of hypothetical interaction in the 
configuration of experiment [20] (see Fig. 1) may be ob- 
tained as the additive sum of the potentials (16) with 
appropriate atomic densities of the lens and the disc ma- 
terials 
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Here p' i: V( (i — 1,2,3) are the densities and volumes 
of the lens and the covering metallic layers {pj, Vj are 
the same for the disc). The proton mass m p appeared 
due to the use of usual densities instead of the atomic 
ones. In numerical calculations of Sec. V the values p[ = 
2.23 g/cm 3 , Pl = 2.4g/cm 3 , p' 2 = p 2 = 8.96 g/cm 3 , p' 3 = 
p 3 = 19.32 g/cm 3 [20] will be used. 

The hypothetical force between lens and disc can be 
computed as the derivative 
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The integration (17) can be performed most simply in 
a spherical coordinate system described in Sec. II (see 
Fig. 1). In these coordinates the multiple integral from 
(17) takes the form 
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The quantities 6\, # 2 were defined in (4), the other inte- 
gration limits are (as above by the prime the lens param- 
eters are notated) 
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For generality the thicknesses of metallic layers on the 
lens and on the disc are permitted to be different. 

To calculate the integrals in (19) it is convenient to use 
the expansion into a series of spherical harmonics [34] 
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where I v {z), K v (z) are Bessel functions of imaginary ar- 
gument. For large arguments (z » 1) one has 
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Let us consider separately the cases of small and large 
parameter A. In (21) the condition z ^> 1 corresponds to 
A -C r\,T2 which is valid when A <C R (actually A must 
be less than the lowest size of the interacting bodies, i.e. 
A < H in the configuration under consideration). Sub- 
stituting the asymptotics (22) into the right-hand side of 
(21) and taking into account the completeness relation 
for spherical harmonics (A9) we obtain for A < H : 
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Substituting (23) into (19) we get two different situ- 
ations depending on the value of index i. If i = 2,3 
(the integration is over the layers covering the lens) the 
integration limits in n do not depend on d: 
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where 9 = min (9\, 6* 2 ). 

It is convenient to consider the force fy'u ( a ) = 
—dUy^ (a) I da instead of the potential energy (24) . This 
gives the possibility to remove the integration with re- 
spect to r 2 . Integrating with respect to the variables r\ 
and i? in a standard way [35] one gets for i = 2,3 and 
.7 = 1,2,3: 

&u (a) = 47r 2 A 3 e- a -/ A (i?- A) (l - e" A *-i/ A ) (25) 
x |l - e - BR/x - e ~ d ' /x 1 - e -(R+ d j) B / x j . 

Here the following notations are used: 
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Now we consider the second situation where the inte- 
gration in (19) is over the lens itself (i = 1). Here the 
quantity Uy^ (a) has the same form as in (24) but with 
integration limit R[ min depending on Considering the 
force instead of potential energy and performing the in- 
tegration by the use of tabulated formulas [35] we find 
the result: 
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f£j\a) =4^ 2 A 3 e -^/ A {( J R-A) 



(27) 



1 - er BR l x - e 



-(R+dj)B/\^ 



-H/X 



A 1 -A 2 e- BH / x -A' 3 -e 



(o) P -dj/X 

3 e 

+A U) e -d J /X e -(H+d ] )B/x'j | _ 

Here the following notations are introduced: 

an = a + A 1 + A 2 + A' x + A' 2 , 
ai2 = a + A 2 + A; + A 2 , 
da = a + A' x + A' 2 , 



(28) 



R 
R 



(R-H)(H + X) 



A 2 = - 
A 



H 

(R-H)(H + Xt Q ) 



H 

j ' R + d j 



X 



R + dj 
H + dj 



H 

(R — H)(H + dj + A) 
(R-H){H + d 3 + Xt 



-A 



A 



to is defined in explanations to (6). 

The complete value of the hypothetical force, according 
to (17), (18), may be presented in the form: 



he 



= - a ^L PiPif: 



Yu ! 



(29) 



where are defined in (25), (27). The expressions 

(25), (27), (29) will be used in Sees. V,VI for calculating 
constraints on Yukawa-type hypothetical interactions. 

Note that for the extremely small A ~ a the expression 
(29) may be additionally simplified: 



F Yv = _ a ^ 4 7r 2 A 3 e- a / A i? 



mi 



(30) 



p[e 



3 -(A' 1+ A' 2 )/X + f -A' 2 /X 



P' 3 



pie -(A 1+ A 2 )/X +p2e -A 2 /X 



■ P3 



Exactly this result (for equal thicknesses of the covering 
metallic layers) follows in the limit of small A from the 
formula derived in [21] for the configuration of a lens and 
an infinite disc. Thus for A ~ a the finiteness of the disc 
size does not influence the value of the hypothetical force. 
At the same time for larger A it is necessary to take into 
account the finite sizes of the disc (unlike the case when 
we calculated the Casimir force in Sec. II). 

Let us start with the case of large A (A 3> R). Now 
we may neglect the exponent in Eq. (16). As a result 
the interatomic potential does not depend on A any more 
and behaves as l/r\2- For this potential the hypothetical 
force between a lens and a disc was calculated in [9] un- 
der the assumption that the disc area is infinitely large. 



Such suggestion is not reliable for the potential under 
consideration due to its slow decrease with the distance. 
The general method developed in the Appendix for the 
potentials of the form r 12 ' 2fc+1 ' ) with k > 1 also should 
be modified to include the case k = 0. We will take into 
account that for large A the covering metallic layers prac- 
tically do not contribute to the result. By this reasoning 
one may integrate directly over the lens and the disc vol- 
umes. As a result in a spherical coordinate system used 
above the expression for the hypothetical force is: 

F Yu~~a^p'\P\ J dip i fdip 2 Jsiwdxd-dx (31) 
p o o o 

02 R -Rmax 

x J sin d 2 dd 2 J r\ dr x j r\ dr 2 



r 2 cos -$2 — Tl cos$i 



The calculational details for integration in (31) are 
given in Appendix. The result is: 



tic 

F Yu w -a^p[ Pl 2ir 2 L 2 DH 



L 2 ( D + H H. D + H 
1 ttttt I m 1 In 



(32) 



4RH 



D D 



H 



In the intermediate range between A < H and A 3> R 
the integration in (17) was performed numerically For 
this purpose algorithm 698 from netlib [36] was used. It 
is an adaptive multidimensional integration, the FOR- 
TRAN program is called DCUHRE. A large number of 
function calls (about 5 x 10 6 ) for A > 2.3 x 10~ 4 m was 
necessary in order to obtain reliable results. For smaller 
values A the program does not work. The results of nu- 
merical calculations are in good agreement with the an- 
alytical ones (see Sees. V,VI). 



IV. CALCULATION OF A POWER-LAW 
INTERACTION 

In this Section the force is calculated which may act 
in the configuration of experiment [20] due to power-law 
hypothetical interactions. The power-law potential be- 
tween two atoms belonging to a lens and a disc is 



V n = -X n N l N 2 Uc 



±(n L \ 



r\ 2 \ r i2/ 



(33) 



Here A„ is a dimensionless constant and r = IF = 
10~ 15 m is introduced for the proper dimensionality of 
potentials with different n [18]. 

The potential energy of the lens and the disc (see 
Fig. 1) due to a hypothetical interaction is an additive 
sum of the potentials (33) with appropriate atomic den- 
sities. It can be written in a form analogous to (17) 
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h ro— l 3 

p ij=l 



(34) 



F n (a) = - 



dUn(a) 
da 



The quantities Un here are dehned by the same inte- 
grals as in (17) where instead of cxp(— ri 2 /A) the function 
l/r™2 _1 should be substituted. 

At first we consider the case n — 3. Here the result 
for ujf'^ is given by (All) with k = 1. We rewrite it 
by using (A6) for r\r 2 and introducing the new variables 
t = cosi9, x\ = tri, x 2 = tr 2 



1 dt %r dxTr*K {Xl+X2) \ (35) 



x 2 - xi 



tR'. 



where to is defined in the explanation to (6). 

Calculating the force, one can eliminate the integration 
with respect to x 2 : 



f j \a) = 



da 



tR', 



"7? I 



dxi 



(xx + tR 



j, max 



(36) 

) 2 (X! +ii? J , min ) 2 " 



tR 



to tR'. 



j, max 



Xi 



tRj : 



Xi 



Performing the integration in (36) for different val- 
ues of i, j to lowest order in the small parameters H/R, 
L 2 /(2R 2 ), D/R, A'JH, Aj/D, a/H, a/D and a/L, and 
substituting the results into (34), we obtain: 



F 3 ( a ) = -47r 2 A 3 ^i? 



mi 



(37) 



x \ PiPi 



+ Pi 



, 2RD + L 2 L 2 , (2RD + L 2 )H 

D In 1 In — = — 

2RD 2R L 2 {H + D) 

L 2 



p 2 Ai In 



+/9 3 A 2 In 
Pi ( P2K In 



2R(a + A 1 + A 2 + A; + A' 2 ) 
L 2 



2i?(a + A 2 + A' 1 + A' 2 ) 
D 



a + Ai + A 2 + A' 2 



+p' 3 A' 2 In 
+ p 2 ( ^A; In 



D 



a + Ai + A 2 
a + Ai + A 2 + A 2 



+p 3 A 2 In 



a + A 2 + A 2 

a + Ai + A 2 



a + A 2 



+p 3 p 2 A'j In 



a + A2 A +A2 +^A-ln: 
a + A; 13 2 



Here L 2 = min(L 2 , 2RH). Note that in the specific case 
Ai = A 2 = Ai = A' 2 = A, p' 2 = p 2 , p' 3 = p 3 , and L = L 
Eq. (37) coincides with Eq. (16) of Ref. [22]. 

Let us discuss now the power-law potentials with the 
even powers n — 2, 4. In this case the corresponding 
quantities U 2 ^ (k — 1,2) from Eq. (34) are most con- 
veniently represented in cylindrical coordinates (r, ip, z) 
with the origin at the lens top and the z-axis orthogonal 
to the disc surface. Then the quantities U^J^ describing 
the interaction energy of the lens atoms with the atoms 
of the disc and its covering layers take the form 

H /Oi) —aj L 

U^ k ' j) (a) =4tt J dz x J dn J dz 2 j dr 2 (38) 

A'j+A^ -aj-dj 
Ti T 2 dip 



I 



[t 2 + t| - 2riT 2 cos <p+(zi- z 2 ) 2 ] 



In (38) the following notations are introduced 

f(z 1 )=[2(R-A' 1 -A' 2 )z 1 -z 2 ]-\ (39) 
ai=a + Ai+A 2 , a 2 =a + A 2 , a 3 = a, 

the thicknesses dj were defined in (26). 

(i j) 

The quantities U 2k with z = 2,3 describe the interac- 
tion energy of the lens covering layers with the disc and 
its layers. They are expressed by 

U^\a) =4tt J sintfidtfi J dn J dz 2 J dr 2 (40) 

R'. . -ai—di 



7T 

/ 



r 2 r 2 dtp 



t 2 + r| - 2tit 2 cos <p + (21 - z 2 ) 2 ] 



Here the layers covering the lens are described in spher- 
ical coordinates used above, so that n = rising, 
z\ = R — n cosi?i. 

To perform the integration with respect to <p in (38), 
(40) it is helpful to use the integral representation of Leg- 
endre polynomials [35] . Also the integration with respect 
to z 2 can be removed when calculating the force instead 
of the potential energy. The result for i = 1 is 

H f(zi) L 

F 2 [ l' i \a)= J dz x fr 1 dnfr 2 dT2 (41) 





A' 1+ A' 2 



In the same way for i = 2, 3 the result is 
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7 2k\ a ) = j Sinl? l M l j r l dT l j T 2 dT 2 



sini?i d$i 

R 

2 



(42) 



x 



5=i (4?) fc 



In (41), (42) the following notations are introduced 



» 



(p)\2 



^+r^ + (z 1+ ^) 



4t 1 t 2 j 5 ) 



1 



(p)\2 



r/+^ + (z 1+ ^) 



(43) 



The quantities with i = 2, 3 are obtained from 

, by substitution of n , zi according to explana- 
tions after Eq. (40). 

Eqs. (41), (42) with the notations (43) look rather cum- 
bersome. In spite of this all involved integrals can be cal- 
culated explicitly by the use of Ref. [35] . As a result the 
hypothetical force Fi(a) computed according (34), (41), 
(42) in the lowest order in small parameters mentioned 
above is: 



F 4 (a) 



-2tt 2 A 4 tt-R 



mt 



(44) 



PiPi 



In- 



2R(a + Ai + A 2 + A; + A' 2 ) 2RH 
, ( , a + Ai + A 2 + A; + A' 2 

+ P 1 \ P2 In t T-; T-. 

1 V a + A 2 + Ai + A 2 

a + A 2 + Ai + A 2 



+P3 In ■ 



+ Pi ( P'i 

+ P2 



a + Ai + A 2 
Ai 



a + Ai + A 2 + A' 2 ' Ps a + A , f A 2 
Ai ( Ax 



A 2 



a + A 2 + A 2 V P2 a + Ai + A 2 + A 2 
A 2 



, , A 2 
+P 3 " 



+P3- 



P2 



Ai 



A, 



a + A 2 V a + Ai + A 2 



P3- 



The hypothetical force F 2 (a) calculated by (34), (41), 
(42) in the lowest order of the same small parameters is: 



^-D>Rln 2RD + L2 



-^{L 2 R-D 2 R-DL 2 )\n 
+ L^ln^ + L^ln 2 ^ + ? (45) 



2RD 

2RD + L 2 
I 5 



H 



2DR + L 2 
L 4 



+ \w {L2 + DR)ln mH-2-nn-m 



Note, that here only the contribution of the lens and 
the disc materials is written out. The covering metallic 
layers practically do not contribute to the value of the 
force for the power-law potentials with n = 2 [22] . 



V. CONSTRAINTS FOR HYPOTHETICAL 
INTERACTIONS FROM THE RECENT 
EXPERIMENT 

The results of Sees. II-IV are used here for obtain- 
ing stronger constraints on the constants of hypothetical 
Yukawa- and power-type interactions which follow from 
the measurements of experiment [20] . The absolute error 
of the force measurements in [20] was AF — 1 ^idyn = 
10~ n N in a range of distance between the lens and the 
disc from a = 1 /j,m till a = 6 /zm. With this error the 
expression (1) for the Casimir force was confirmed and 
no corrections to it or unexpected interactions were ob- 
served. 

We now discuss the values of a for which the most 
strong and reliable constraints on hypothetical interac- 
tions can be obtained from the above mentioned result. 
Thereby the corrections to the Casimir force considered 
in Sec. II have to be taken into account. Evidently, be- 
cause all kinds of hypothetical forces decrease with dis- 
tance the smallest values of a are to be prefered. From 
this point of view, e.g., a = 1 ^m should be chosen. But 
for these values of a, as it follows from Sec. II, the theo- 
retical value of the force under measuring is not strictly 
defined. Although the Casimir force itself is rather large 

(F^ w —3.1 x 10~ 10 N) and the temperature correc- 



tion to it is rather small (AtF ( 



(o) 



-0.8 x 10~ n N) the 



corrections due to finite conductivity and due to surface 
distortions are large. According to Fig. 2 (curve 1) the 
corrections to finite conductivity As Fq^ at a = 1 /im 
is of order -0.23F, 



(o) 
c 



7.1 x 10~ n N. The correction 
due to short-scale distortions A^F^ (curve 2) is of 
order 0.12F^ 0) w -3.7 x 10 _11 N. At the same time 
the correction due to large-scale deviations of bound- 
ary surfaces from the perfect shape AiF^ may achieve 

±0.3F£ 0) w t9.3 x 10" n N. All these corrections are of 
the same order or much larger than the absolute error of 
force measurements AF. Moreover, the largest correc- 
tion due to large-scale deviations can not be estimated 
theoretically because the actual shape of interacting bod- 
ies was not investigated in experiment [20] . In this situa- 
tion the cancellation of contributions from different cor- 
rections to the force value occurs very likely. 

The constraints on the parameters of hypothetical in- 
teractions a, A, A„ may be calculated from the inequality 



Fp{a)+F h {a)-Fg>{a)\<AF, 



(0), 



(46) 



where Fq 1 is the theoretical Casimir force value with ac- 
count of all corrections, Fh is the hypothetical force Fy u 
or F„ calculated in Sees. Ill, IV. 
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Substituting the general expression for Fq 1 into (46) 
one obtaines 

\F h (a) + A E (a) + A, 4%)! < AF, (47) 



where 

A E (a) = A T 4 0) (a) + A, F^ 0) (a) + A d F^ (a) (48) 

is the sum of corrections whose values for different values 
of a are accessible. 

Although the values of the correction due to large-scale 
deviations are unknown its dependence on a is given by 
Eq. (12). Thus for two different values of a it holds 

A*4% 2 ) = ±A lF ^( ai ) 7 k 21 = g)\ (49) 

According to the results of Sec. Ill the value of the 
hypothetical force is proportional to interaction constant 
Fh(ai) = aG\(di), or to A„ for power-law interactions 
(see Sec. IV), with some known functions G\ (G n ). Con- 
sidering (47) for two different values of distance 0.1,0.2 
with account of (49) and excluding the unknown quan- 
tity AiFg\ai) we obtain the desired constraints for a 
(An) 

-(fei + 1)AF - A E (oi) + A: 21 A s (a 2 ) 

<a[G A (ai)-fc 2 iG A (a 2 )] (50) 
< (fc 2 i + 1)AF - A E (oi) + fc 2 iA s (o 2 ). 

The specific values of a\ , o 2 in (50) should be chosen 
in the interval 1 /im < a < 6 /im in order to obtain the 
strongest constraints on a, X n . For the upper limit of 
the distance interval (a w 6 /im) the Casimir force Fc (a) 
from (14), i.e. together with the temperature correction, 
should be considered as a force under measuring. All 
corrections to it are much smaller than AF = 10~ n N. 
By this reason for such values a the constraints on the 
hypothetical interaction may be obtained, instead of (47), 
from a simplified inequality: 

\F h (a)\ = \aG x (a)\ < AF. (51) 

Now let us turn to numerical calculations of constraints 
starting with the Yukawa-type potential. The constraints 
on a were obtained from Eq. (50). For every A some 
pair of distances ai,a 2 was selected which provides us 
with the strongest constraints. The hypothetical force 
aG\{ai) was calculated by Eqs. (29) for small A, (32) for 
large A and numerically in the intermediate region ac- 
cording to Sec. III. The results are presented in Fig. 3 
by the solid curve 1 which corresponds to attractive hy- 
pothetical force (a > 0) and by the dotted curve 2 cor- 
responding to repulsion (a < 0). In the (A, a)-plane the 
region above the curves is prohibited and the region be- 
low the curves is permitted. Let us consider firstly the 
case a > 0. Here the constraints for A ~ 10~ 7 m were ob- 
tained by using Oi = 1 /im, a 2 = 1.5 /im. With increase of 
A the best constraints follow from a\ = 1 /im, a 2 = 2 /im. 
For A > 8 x 10~ 6 m the values a\ = 1 /im, a 2 = 3 /im were 
used. Now let us turn to the case a < 0. Considering 
A w 10~ 7 m, we have chosen a\ — 1 /im, a 2 = 1.5 /im 
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once more. Then in the range till A w m the val- 
ues a 2 = 2 pm or a 2 = 1.5 pm were used (ai = 1/im). 
For A > 10~ 5 m the system of inequalities (50) does not 
lead to better results than the single inequality (51) used 
with a — 6/xm. The complicated character of curves 1,2 
in Fig. 3 (non-monotonic behaviour of their first deriva- 
tives) is explained by the flaky structure of test bodies. 
For A > 10~ 5 m the metallic layers do not contribute es- 
sentially to the value of the force and it is determined 
mostly by quartz. But for A ~ 10~ 5 m the contribution of 
the layers becomes the main one. 

In Fig. 3 the known constraints are shown also fol- 
lowing from the former Casimir force measurements be- 
tween dielectrics [4-7] (curve 3), Cavendish-type exper- 
iments (curve 4 [16], curve 5 [14], curve 6, a > [14], 
curve 7, a < [15]) and Eotvos-type experiment [17] 
(curve 8). It is seen that the new constraints follow- 
ing from [20] are the best ones within a wide range 
2.2 x 10~ 7 m< A < 1.6 x 10" 4 m. They surpass the old 
ones up to a factor of 30. For A < 2.2 x 10~ 7 m the old 
Casimir force measurements lead to better constraints 
than the new one. This is caused by the smallness of 
the Casimir force between dielectrics comparing the case 
of metals and also by the fact that there the force was 
measured for smaller values of a. 

Now let us obtain constraints for the constants A n of 
power-law hypothetical interactions which follow from 
the experiment [20]. Here the inequalities (50) should 
be used once more, where a has to be replaced by 
A n and G\(cii) by G n (a,i). For the power-law interac- 
tions the hypothetical force A n G„(ai) is calculated by 
Eqs. (37), (44), (45) of Sec. IV. The strongest constraints 
are obtained for a\ — 1 pm, a 2 — 3 /im. Thus for the po- 
tential (33) with n — 2, 3, 4 the new constraints are: 
A 2 < 1.1 x 10- 26 , A 3 < 1.6 x 10~ 14 , A 4 < 3.6 x 10" 3 . 
Note that these constraints are not so strong as the 
ones obtained from the old Casimir force measurements 
(A 2 < 5 x 10" 28 , A 3 < 5 x 10" 15 , A 4 < 3 x 10~ 3 [9]) or as 
the best ones obtained up to day from the Cavendish-type 
experiment (A 2 < 7xl0" 30 , A 3 < 7xl0~ 17 , A 4 < lxl0~ 3 
[19]). The reasons for this are the same as for Yukawa- 
type interactions with small A. Actually, the power-law 
potential with n = 2 leads to the force (45) between a 
lens and a disc which does not depend on a. The depen- 
dence of F3 on a is also very week (see (37)). By this 
reason the covering metallic layers do not contribute to 
the value of the hypothetical force and can not compen- 
sate the factors mentioned above. For n = 4 there is the 
noticeable dependence on a in (44). By this reason the 
constraints obtained from the new experiment is almost 
of the same strength as from the old one. 

As it is shown in the next Section the new constraints 
obtained from the experiment [20] can be considerably 
strengthened owing to some modification of its parame- 
ters. 



VI. POSSIBLE IMPROVEMENT OF THE 
OBTAINED RESULTS 

The experiment |2(i| ] was designed to demonstrate the 
Casimir force between metallic surfaces. The strengthen- 
ing of the constraints on the constants of hypothetical in- 
teractions was derived from its results afterwards. There- 
fore it is likely that modifications of the design would al- 
low to get much stronger constraints. The simplest sug- 
gestion follows from the fact that the hypothetical forces 
are proportional to the product of the densities of the lens 
and the disc (see Sees. Ill, IV). The value of the Casimir 
force measured in [^0| was determined by the thin cov- 
ering gold layers only. For the hypothetical forces which 
decrease with distance not so fast this is not the case. 
By this reason the contribution of the hypothetical forces 
may be increased by the use of some high density metals 
as material for the lens and the disc instead of quartz. 
As such material iridium (p = 22.4 g/cm 3 ) looks very 
promising. The gold cover of 0.5 pm thickness should 
be preserved due to its good conductivity. With a lens 
and a disc made of iridium the values of the hypothetical 
forces (if any) are increased approximately by a factor of 
10 2 . An increase of the lens and the disc geometrical pa- 
rameters (R, H, L, D) to become larger than the values 
used in the experiment p0| ] is not required. According to 
our estimates this would not lead to an essential further 
strengthening of the constraints (because the additional 
volume is situated too far from the nearest points of the 
bodies). 

In Fig. 3 (curve 9 for a < 0, curve 10 for a > 0) 
the constraints for Yukawa-type interactions are shown 
which would follow from an experiment like |20| ] with the 
lens and the disc made out of iridium covered by a gold 
layer. The constraints were obtained using the inequality 
( p0| ) in the same way as in Sec. V for quartz test bodies. 
The hypothetical force was calculated according to Sec. 
Ill with pi — p[ = p. It is seen that the prospective 
constraints of curves 9,10 are about 100 times stronger 
than that obtained actually (curves 1,2) in the range A > 
10" 3 m. For smaller A the strengthening is not so high 
because the gold covering layers contributed hypothetical 
forces essentially in this range in f^o). For A510~ 6 m 
the prospective constraints are approximately the same 
as the actually obtained from [g0| (here the gold layers 
themselves determine the result). It is noticeable that 
with iridium test bodies the Casimir effect would give the 
best constraints for the wider A-range 2.2 x 10~ 7 m < A < 
3.2 x 10~ 3 m and would exceed the constraints following 
from the Cavendish-type experiment of Ref. (curve 
4). 

Better constraints may also be obtained on the con- 
stants of power-law hypothetical forces by use of irid- 
ium test bodies. Performing the calculations using the 
inequalities (|5^) with a\ = 1 /im, a 2 = 3 pm and Eqs. 
([37|), ( [i"5l) for the hypothetical forces one obtains A 2 < 
1.2 x 10 -28 , A 3 < 1.8 x 10~ 16 . In the same way for n = 4 
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with ai — 1 fim, a-i = 2.5 fim one has A4 < 8.4 x 10~ 5 . All 
these constraints are stronger than those obtained before 
from the Casimir force measurements between dielectrics 
(see Sec. V). For n — 4 the prospective constraint is 
stronger than the best one obtained up today from the 
Cavendish-type experiment (A4 < 1 x 10~ 3 ). 

Now let us consider one more modification of the ex- 
periment |2(i(| which causes a further strengthening of the 
constraints. The experiment |2(]] was not so sensitive as 
it might be. This was because of the missing vibration 
isolation from the surrounding building. Using an ap- 
propriate isolation and larger distances between the test 
bodies the absolute sensitivity might reach the 0.01 mi- 
crodyne level |37|] . 

Here we estimate the prospective constraints which 
might follow from the Casimir force measurements be- 
tween the iridium lens and disc covered by gold layers, 
assuming AF = 10" 13 N. We consider the upper limit of 
the a-interval, a = 6 fim, and suggest that the theoret- 
ical value of the Casimir force (U4) including the tem- 
perature part is confirmed with the absolute error AF. 



(Note that for a = 6 /im we have F ( 
and A T F^ } 



(o) 
c 



-1.43 x 10" 12 N 



-2.48 x 10 _12 N, so that the temperature 
part can not be considered as a correction any more.) 
The modulus of corrections to the Casimir force due to 
the finite conductivity and surface distortions are less 
than AF at a — 6 /im (the largest of them would be 
A So F c m 0.36 x 10~ 13 N, see Fig. 2). That is why, as the 
first approximation, it is possible to get the prospective 
constraints on the hypothetical force from the inequality 

The calculational results are shown in Fig. 3 by the 
curve 11. It is seen that the new prospective constraints 
overcome almost all results following from the different 
Cavendish-type experiments (except of a small part of 
the curve 7 for a < 0). They may become the best ones 
in a wide A-range 5 x 10 -7 m < A < 5 m. The maxi- 
mal strengthening comparing the results following from 
ipTjf achieves 10 4 times. For the intermediate A-range the 
strengthening achieves several thousand times. For ex- 
tremely small A the promised results are weeker; this is 
caused by the use of inequality ( pT| ) instead of the more 
exact one, Eq. (|5(]). Strictly speaking, the special inves- 
tigation of large-scale deviations of the surface from the 
ideal shape and short-scale distortions is needed to obtain 
the most strong constraints on the constants of hypothet- 
ical interactions from the Casimir force measurements. 

The prospective constraints of curve 11 would give 
the possibility to restrict the masses of the spin-one 
antigraviton (graviphoton) and dilaton. The interac- 
tion constants of Yukawa-type interaction due to ex- 
change by these particles are predicted by the theory: 
a a fn 10~ 40 , a d « 2 x 10" 39 @. As it is seen from 
curve 11 of Fig. 3, e.g., for graviphotons the permit- 



ted values of A are A„ < 4 x 10 



or for its mass 



m a > V(^™ aXc ) ~ 5 x 10 eV. These constraints are 
stronger than those known up to date (A a < 3 x 10~ 3 m, 



m a > 6x 10~ 5 eV |12 ) obtained from Cavendish-type ex- 
periment. Note that obtaining much stronger constraints 
for the parameters of the graviphoton is of special interest 
in connection with the recently claimed experimental ev- 
idence for the existence of this particle from geophysical 
data H§. 

Decreasing of the absolute error of force measurements 
till AF — 10 _13 N will give the possibility to strengthen 
the constraints on power-law interactions as well. For 
a lens and a disc made of iridium and covered by gold 
layers the results are obtained from the inequality ( |5l| ) 
with a = 6 fim: A 2 < 2.13 x 10~ 30 , A 3 < 3.25 x 10^\ 
A4 < 1.8 x 10~ 6 . These constraints overcome essentially 
the current results following from the Cavendish-type ex- 
periment (see Sec. V). The greatest strengthening by 
several hundred times takes place for n = 4. That is 
why the realization of the experiment on demonstration 
of the Casimir force with improved parameters is of great 
interest for the problem of hypothetical interactions. 



VII. CONCLUSION AND DISCUSSION 

In this paper we performed a careful calculation of the 
Casimir and hypothetical forces according to the config- 
uration of experiment [20], i.e. between a spherical lens 
and a disc made of quartz whose surfaces were covered by 
thin layers of copper and gold. The finiteness of the disc 
area was taken into account and the corrections were cal- 
culated to the Casimir force between a lens and a disc of 
the infinite area (Sec. II). These corrections were shown 
to be negligible. That is why the use of the theoretical 
result for infinite disc in [20] for confronting with exper- 
imental data is justified. 

Different corrections to the Casimir force were dis- 
cussed, e.g., due to the finite conductivity of the bound- 
ary surfaces, deviations of their geometry from the per- 
fect shape and due to non-zero temperature (Sec. II). It 
was shown that the corrections due to finite conductivity 
and to short-scale distortions have the opposite sign and 
may partly compensate each other. At the same time 
the global deviations of the boundary surface geometry 
from the perfect shape lead to both positive and negative 
corrections (which may reach 30% of the Casimir force 
acting in a perfect configuration with space separation 
a = 1 fim) . By this reason a detailed investigation of the 
boundary surface geometry is required when confronting 
experimental and theoretical results for configurations 
with a small space separation. As to the temperature 
contribution to the Casimir force it may be considered 
as a correction for the small space separation only and 
should be included into the force under measuring start- 
ing from a w 3 fim. 

The calculation of the Yukawa-type hypothetical force 
has shown that for large values of A~10~ 3 m it practi- 
cally does not depend on a, for A ~ 10~ 5 m the covering 
metallic layers do not contribute to its value essentially, 
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but for smaller A the contribution of the layers becomes 
the main one (Sec. III). For the power-law interactions 
the strongest dependence of the force on a was obtained 
for n = 4 (Sec. IV). For n = 2 it is practically absent and 
for n = 3 there is only a weak dependence of the force on 
a. As a result the covering metallic layers practically do 
not contribute the value of force for n = 2, 3. 

The careful calculation of the Casimir and hypotheti- 
cal forces in the configuration of experiment [20] gave the 
possibility to obtain reliable constraints on the param- 
eters of hypothetical interactions (Sec. V). In the case 
of Yukawa-type hypothetical interactions the new con- 
straints surpass the old ones following from the Casimir 
force measurements between dielectrics by a factor of 30 
in a wide range 2.2 x 10~ 7 m < A < 1.6 x 10~ 4 m (curves 
1,2 in Fig. 3). In this A-range the obtained constraints 
are the best ones on the Yukawa-type interactions fol- 
lowing from laboratory experiments. For the power-law 
interactions the experiment [20] does not lead to new con- 
straints which would be better than the ones known up 
to date. 

According to the analysis presented above, by some 
modification of parameters of the experiment [20] the re- 
lated constraints on the constants of hypothetical inter- 
actions may be essentially improved (Sec. VI). With the 
use of iridium test bodies (instead of quartz ones) the 
possible improvement is shown by the curves 9, 10 in 
Fig. 3 and achieves two orders of magnitude. If in ad- 
dition the accuracy will be improved due to the use of 
vibrational isolation from the surrounding building the 
resulting constraints are given by the curve 11 in Fig. 3. 
In this case the improvement will achieve four orders of 
magnitude for large A and several thousand times for in- 
termediate A. It is notable that the improved Casimir 
force measurement promises more strong constraints for 
the range 1.6x 10~ 4 m < A < 5 m than the ones known up 
to date from the Cavendish- type experiments. Obtaining 
such constraints would supply us with new information 
about light hypothetical particles, e.g., graviphoton, dila- 
ton, scalar axion etc. The experiment with iridium test 
bodies and improved accuracy will give the possibility 
to strengthen constraints on power-law interactions as 
well. Comparing the current constraints following from 
the Cavendish-type experiment the prospective strength- 
ening is by a factor of 3.3, 21.5 and 555 for n = 2,3,4, 
respectively. 

To conclude, we would like to emphasize that the new 
measurements of the Casimir force are interesting not 
only as the confirmation of one of the most interesting 
predictions of quantum field theory (there are also other 
experiments on Casimir effect in preparation, see, e.g., 
[39]). The additional interest arises from the hope that 
it may well become a new method for the search of hypo- 
thetical forces and associated light and massless elemen- 
tary particles. 
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APPENDIX 

In this Appendix we present the derivation of sev- 
eral mathematical expressions used in Sees. II-IV. Let 
us start with the calculation of the multiple integrals 



27T 0! 



hk+i = J dtp 1 Jsiwd 1 d-d 1 j r\dr x (Al) 

fl»i„(*i) 

x j dtp 2 j sin i? 2 dd 2 J r\dr 2 r^ 2 2k ^ x 



o o iJ mi „(# 2 ) 

whose integration limits satisfy the conditions 



R min 



mm ^ -"-max: 



sin 0i < 1, sin0 2 <Cl, 



-Rmax -^min ^ -^min' 



(A2) 



Integrals of that type were essential in Sec. II (for k = 3) 
and in Sec. IV (for k = 1). 

To calculate I 2 k+i it is convenient to use the following 
expansion into a series of spherical harmonics [34] 

rr-f- 1 = 4^]T -^- a j- (2fc+1)1 (ri,r 2 ) (A3) 
/ 

X £ ^(01,¥>l)*W<>2,¥>2), 



where the radial part may be represented in the form 



[-(2fc+l)]/ \ 



(|), (ri+r 2 ) 2l + 2k+1 



(A4) 



xf(l + ^l il+1;2i + 2; ^_ 

2 {n+r 2 y 



Here F(a, (3; 7; z) is the hypergeometric function, (n); = 
T(l + n)/T(n), T(z) is the gamma function, and in accor- 
dance with inequalities (A2) r 2 > r\. 
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It is readily seen that due to (A2) the argument z = 
4rir 2 /(ri + r 2 ) 2 of the hypergeometric function in (A4) 
is of order of unity. So according to [35] it is profitable 
to use the hypergeometric function of the argument z\ — 
1-z: 

2k + 1 \ r(2/ + 2)r(±^) 

F[l + —^—,l + l;2l + 2;z^- V ' V 2 ; 



3-2k 



2k + 1 2k + 1 

x F{ 1 + — - — ,1 + 1; — - — ;zi 



+ 



T2 + n 



2 ' ' 2 
2fe - x T(2l + 2)T( 2 ^) 
T(l+ 2Mi) T{1 + 1) 



(A5) 



^ . , 3 - 2k , 1 3 - 2fc 
x F| / + — - — ,Z + 1; — - — ;zi 



Note that due to inequalities (A2) it holds « 1. On 
account of this the first contribution on the right-hand 

fe-i 

side of (A5) is of order z x 2 relatively the second one and 
it is possible to omit it for k > 1 . One can also substitute 
the hypergeometric function of the small argument z\ in 
the second contribution to (A5) by unity. In addition the 
product of radiuses in (A4) may be expressed in terms of 
their sum with the same accuracy 

ri r 2 = 1 [(n + r 2 f - (n - r 2 ) 2 ] w 1 (n + r 2 ) 2 . (A6) 

Substituting (A5) and (A6) into (A4) and using the 
properties of gamma function we obtain 



[-(2fc+l)], x 



2(2/ + 1) 



(2k- l)(n +r 2 ) 2 fa-n) 2 *- 1 ' 

(A7) 

As a result expansion (A3) with the condition k > 1 
takes the form 



-2fc-l 



87r (r 2 - n 



,l-2fc 
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2fc - 1 (n + r 2 ) 2 

oo I 

X E E ^(^l)W^ 2 ). 
;=0 m=-J 



(A8) 



After the substitution of (A8) with the use of com- 
pleteness relation for the spherical harmonics [34] 



oo / 



;=0 m=-l 

x <5(cos^i - cosi? 2 ) (A9) 
Eq. (Al) may be represented as 

2n 2ir 0i 



hk+l 



8tt 



2k -1 



o o 



i 5{ipx - tp 2 ) J sin??i d$i 



< / sini9 2 rfi? 2 <5(cosi?i — cosi? 2 ) 
o 

/ dri / cir 2 



(A10) 



/ / 



r 2 r 2 
'1 '2 



(n + r 2 ) 2 (r 2 - n) 



_ r , \2k-\ ' 



After the integration with 5-functions in (A10) the re- 
sult is 



hk+l — 



16tt 2 
2k- 1 



min(#i ,#2) 



sin t? dti 



(All) 







x J dri J dr 2 

W . (1?) Rmin(#) 



'1 '2 



(ri + r 2 ) 2 (r 2 - n) 



,2fc-l ' 



Eq. (All) is useful for the calculation of the Casimir 
force (Sec. II) and power-law hypothetical interaction de- 
creasing as the third power of distance (Sec. IV). 

Now let us calculate the integral (31) which expresses 
the asymptotic behaviour of Yukawa-type interaction for 
large A (Sec. III). We use once more the expansion (A3) 
into the spherical harmonics with k = 1. For the radial 
part, instead of (A4), it is more convenient to apply the 
equivalent representation [34] 



(-3)/ \ 



21 + 1 r[ 21 + 1 



l-rlr l 2 +1 r 2 (r 2 2 -r\) \r 2 



r 2 



n 



(A12) 

To lowest order in the small parameters H/R, D/R 
it holds ri/r 2 cos^ 2 . Note also that only the term 
with m = from (A3) gives non-zero contribution when 
integrating with respect of <pi. 2 in (31). As a result (31) 
may be rewritten as 

Tic f 

F Yu ~ -a — pi p x I sint? 2 d-& 2 <I>(cosz9 2 ), (A13) 

m p J 



where the following notation is introduced 

0i R 
$(cosi? 2 ) = 47r 2 J sin^i ddi j r\dr\ 



R+a+D 

cos 



/ 



r 2 dr 2 



r 2 cosi? 2 — n cos$i 



a+q 

cos tf 2 



x ]T(cos tf 2 ) ; (2/ + 1) Pi (cos tfi) P ; (cos i? 2 ), (A14) 
and Pi(-z) are Legendre polynomials. 
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Now it is useful to change variables in (A13), (A14) 
according to t\ = cosi?i, t 2 = cosi9 2 , x 2 — t 2 r 2 . After 
that (A13), (A14) take the form 



I „ : -« ^| Pi Pi j dt 2 $(t 2 ), 



mt 



(A15) 



COS $2 



where cosfl 2 = R/VR 2 + L 2 w 1 - L 2 /{2R 2 ) and 



$(£ 2 )=47r 2 y dti jrjdri J x 2 dx 



R+a 



x 2 - n ti 

2 2 _ 2 j.2 
x 2 ' 1 t 2 



^4(2/ + l)fl(ii)^(i2), 



(A16) 



z=o 



with cos6»i = (R-H)/R. 

It is not difficult to calculate the integral (A15) approx- 
imately taking into account that L 2 /(2R 2 ) <C 1 and ex- 
panding <&(t 2 ) into a Taylor series near the point t 2 = 1: 



F^w-a^^i J[$(l) + $(cos0 2 )] 1 
m 2 2 



-a p 1 p x 



1 7 2 
$(1) - -$'(1) — 5 
W 2 W 2i? 2 



2i? 2 
2i? 2: 



(A17) 



where $'(1) = (d*(t 2 )/dt2)|t 2 =i. 

The value $(1) can be calculated with account of 
equality 

00 

£(2i + l)fl(ti)fl(t 2 ) = 2<5(ti-t 2 ) ) (A18) 

1=0 

which follows from the completeness relation (A9). The 
result for the lowest order in small parameters H/R, D/R 
is 



$(1) =4tt 2 R 2 DH 



(A19) 



Now let us find the contribution of the second term 
from the right-hand side of (A17). For this purpose we 
differentiate (A16) 



R+a+D 



$'(1) = 4tt 2 J dtx j r\ dn J x 2 



dx 2 



(A20) 



COS#i R—H 



R+a 



2(x 2 -rih)rl 



K x 2 ~ r l) , n 



+ 



x 2 - n *i 



x 9 r-i 



— £4 (2/ + l)P ; (ii)^fe) 



i=0 



The first contribution on the right-hand side of (A20) 
can be calculated by the use of the completeness relation 



(A18) and for the lowest order in H/R and D/R this 
results in 

It is easily seen that in the lowest order in H/R and 
D/R the second contribution to (A20) is proportional to 
2tt 2 H 2 DR, i.e. is a quantity of the order H 2 /R 2 com- 
paring the first one (A21). By this reason the second 
contribution on the right-hand side of (A20) may be ne- 
glected. 

Substituting (A19) and (A21) into (A17) we come 
to the result (32) for the asymptotic behaviour of the 
Yukawa-type hypothetical interaction in the limit of large 
value of A. 



[1] H.B.G. Casimir, Proc. Kon. Nederl. Akad. Wet. 51, 793 
(1948). 

[2] V.M. Mostepanenko and N.N. Trunov, The Casimir 
Effect and Its Applications (Clarendon Press, Oxford, 
1997). 

[3] M.J. Sparnaay, Physica 24, 751 (1958). 

[4] B.V. Derjaguin, I.I. Abrikosova, and E.M. Lifshitz, 

Quart. Rev. 10, 295 (1956). 
[5] D. Tabor and R.H.S. Winterton, Proc. Roy. Soc. Lond. 

A 312, 435 (1969). 
[6] Y.N. Israelachvili and D. Tabor, Proc. Roy. Soc. Lond. 

A 331, 19 (1972). 
[7] S. Hunklinger, H. Geisselmann, and W. Arnold, Rev. Sci. 

Instr. 43, 584 (1972). 
[8] V.A. Kuz'min, I.I. Tkachev, and M.E. Shaposhnikov, 

JETP Letters (USA) 36, 59 (1982). 
[9] V.M. Mostepanenko and I.Yu. Sokolov, Phys. Lett. A 

125, 405 (1987). 
[10] E. Fischbach et al, Metrologia 29, 215 (1992). 
[11] V.M. Mostepanenko and I.Yu. Sokolov, Phys. Lett. A 

132, 313 (1988). 
[12] V.M. Mostepanenko and I.Yu. Sokolov, Phys. Rev. D 47, 

2882 (1993). 

[13] M. Bordag, V.M. Mostepanenko, and I.Yu. Sokolov, 

Phys. Lett. A 187, 35 (1994). 
[14] J.K. Hoskins et al, Phys. Rev. D 32, 3084 (1985). 
[15] Y.T. Chen, A.H. Cook, and A.J.F. Metherell, Proc. Roy. 

Soc. Lond. A 394, 47 (1984). 
[16] V.P. Mitrofanov and O.I. Pononareva, Sov. Phys. JETP 

67, 1963 (1988). 
[17] B.R. Heckel et al, Phys. Rev. Lett. 63, 2705 (1989). 
[18] G. Feinberg and J. Sucher, Phys. Rev. D 20, 1717 (1979). 
[19] V.M. Mostepanenko and I.Yu. Sokolov, Phys. Lett. A 

146, 373 (1990). 
[20] S.K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 
[21] M. Bordag, G.T. Gillies, and V.M. Mostepanenko, Phys. 

Rev. D 56, R6 (1997); D 57, 2024 (1998). 
[22] G.L. Klimchitskaya, E.R. Bezerra de Mello, and V.M. 

Mostepanenko, Phys. Lett. A 236, 280 (1997). 



15 



[23] G.L. Klimchitskaya and Yu.V. Pavlov, Int. J. Mod. Phys. 

A 11, 3723 (1996). 
[24] J. Blocki, J. Randrup, W.J. Swiatecki, and C.F. Tsang, 

Ann. Phys. (N.Y.) 115, 1 (1978). 
[25] V.M. Mostepanenko and I.Yu. Sokolov, Sov. Phys. Dokl. 

(USA) 33, 140 (1988). 
[26] M. Bordag, G.L. Klimchitskaya, and V.M. Mostepa- 
nenko, Int. J. Mod. Phys. A 10, 2661 (1995). 
[27] E.M. Lifshitz and L.P. Pitaevskii, Statistical Physics, 

Part 2, (Pergamon Press, Oxford, 1980). 
[28] V.B. Bezerra, G.L. Klimchitskaya, and C. Romero, Mod. 

Phys. Lett. A 12, 2613 (1997). 
[29] CM. Hargreaves, Proc. Kon. Ned. Akad. Wet. B 68, 231 

(1965). 

[30] J. Schwinger, L.L. DeRaad, and K.A. Milton, Ann. Phys. 
(N.Y.) 115, 1 (1978). 

[31] V.M. Mostepanenko and N.N. Trunov, Sov. J. Nucl. 
Phys. (USA) 42 812 (1985). 

[32] M. Bordag, G.L. Klimchitskaya, and V.M. Mostepa- 
nenko, Phys. Lett. A 20, 95 (1995). 

[33] L.S. Brown and G.J. Maclay, Phys. Rev. 184, 1272 
(1969). 

[34] D.A. Varshalovich, A.N. Moskalev, and V.K. Kherson- 
skii, Quantum Theory of Angular Momentum (World Sci- 
entific, Singapore, 1988). 

[35] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Se- 
ries and Products (Academic, New York, 1980). 

[36] URL: http://www.netlib.org 

[37] S.K. Lamoreaux, Private communication. 

[38] V. Achilli, P. Baldi, G. Casula et al, Nuovo Cim. B 112, 
775 (1997). 

[39] R. Onofrio and G. Carugno, Phys. Lett. A 198, 365 
(1995). 



1G 




FIG. 1. Configuration of a spherical lens and a disc. Point O is the center of curvature of the lens with height H and 
curvature radius R. The thickness of the disc is D and its radius is L, a being the distance between the disc and the lens. The 
thicknesses of Cu and Au layers on the lens are A' 1: A 2 and on the disc — Ai, A2 respectively. The small volume added to the 
disc in calculations is shown by hatching. 
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FIG. 2. The relative role of different corrections to the Casimir force in configuration of a lens and a disc. Curve I shows 
the correction due to the finite conductivity (AggF^'), curve 2 shows the correction due to short-scale distortions (AdF^'), 
curve 3 shows the temperature correction (AtF^). By the curve 4 the Casimir force itself is shown. 
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FIG. 3. Constraints for the constants of hypothetical Yukawa-type interactions following from the force measurements. 
Curve 1 follows from the new measurement of the Casimir force (a > 0), curve 2 shows the same with a < 0, curve 3 results 
from the old Casimir force measurements between dielectric bodies [4] . Constraints from Cavendish- type experiments are shown 
by curve 4 [16], curve 5 [14], curve 6, a > [14], curve 7, a < [15]. Curve 8 follows from the E6tv6s-type experiment [17]. 
Prospective constraints are shown by the curves 9-11: from the Casimir force measurements with iridium test bodies (curve 9, 
a > and curve 10, a < 0), with iridium test bodies and the improved accuracy (curve 11). The regions of (A, a)-plane below 
the curves are permitted and above the curves are prohibited. 
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